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One of the first orders of business for LHC experiments after beam turn-on will be to 
calibrate the detectors using well understood Standard Model (SM) processes such as W 
and Z production and tt production. These familiar SM processes can be used to cali- 
brate the electromagnetic and hadronic calorimeters, and also to calibrate the associated 
missing transverse energy signal. However, the presence of new physics may already affect 
the results coming from these standard benchmark processes. We show that the presence 
of relatively low mass supersymmetry (SUSY) particles may give rise to significant devi- 
ations from SM predictions of Z+jets and VF+jets events for jet multiplicity > 4 or > 5, 
respectively. Furthermore, the presence of low mass SUSY may cause non-standard devi- 
ations to appear in top quark invariant and transverse mass distributions. Thus, effects 
that might be construed as detector mal-performance could in fact be the presence of new 
physics. We advocate several methods to check when new physics might be present within 
SM calibration data. 
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1 Introduction 



The CERN Large Hadron Collider (LHC) is expected soon to begin circulation of beams, with 
proton-proton collisions at center-of-mass energy y/s ~ 10 TeV, followed by a physics run at 
y/s = 14 TeV. This energy scale ought to be sufficient to explore the mechanism behind the 
breakdown of electroweak symmetry (EWSB) PQ. Expectations are high that the upcoming LHC 
era will provide the necessary data to construct a new paradigm for the laws of physics as we 
know them, and new physics beyond that of the Standard Model (SM) is widely anticipated. 
Numerous theoretical constructs have been created, including weak scale supersymmetry[2], 
theories with extra dimensions [3], little Higgs modelsjl], etc. [5]. Which- if any- of these theories 
correctly describes nature at the TeV scale will be determined by LHC experiments [6]. 

However, before new physics searches can commence at the LHC, it will be necessary to fully 
understand the responses of the ATLAS and CMS detectors to particle physics scattering events. 
Towards this end, the experimental groups will first focus their interest on well-understood SM 
processes. These will include the following: 

• Detect familiar di-jet or multi-jet events and compare the measured jet E T , r/(=pseudorapidity) 
and invariant mass distributions to expectations from QCD theory and Monte Carlo event 
generators. 

• Detect the Z — > ££ events {£ = e or \x initially) and verifying that the Z-peak appears 
with the correct mass and width as measured by LEP/Tevatron experiments. One also 
wants to study the associated multi-jet activity and use the rates as a tune for Monte 
Carlo generators. In particular, Z — > ££+jets production should serve as an excellent 
calibrator for Z — > z/P+jets production, which is one of the most important backgrounds 
for SUSY jets +Ef ss searches. 

• Detect W — > £vi events, verify the expected transverse mass distribution and investigate 
again the associated multiplicity of jets expected. 

• Detect tt production and check that the LHC-measured m t value is in accord with the 
most recent Tevatron results. The simplest channel to begin with is t — ^ blvi while the 
other t decays hadronically, thus allowing a direct m t reconstruction with relatively few 
ambiguities. 

Once these familiar processes are measured, and it is ascertained that expected results are in 
accord with previous LEP/Tevatron measurements and SM theory, then the hunt for deviations 
from the SM can be made and the search for new physics will have begun. In addition, if a high 
degree of confidence is gained in the observation of W and Z production, then these reactions 
may serve as a luminosity monitor. 

The rate for boson production at the LHC is ~ 100 nb, while Z production occurs at 
the ~ 40 nb rate. The pp — > tiX cross section for m t = 175 GeV is ~ 800 pb. Thus, with only 
0.01 fb _1 of integrated luminosity, sizable W, Z and tt signals are expected. 

In this paper, we wish to make the point that it is possible that new physics signals at 
the LHC may be so large that they may already make their presence felt in these baseline SM 
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reactions. In this case, they would potentially disrupt efforts at calibrating both detectors and 
Monte Carlo event generators to expected SM signal levels. 

As an example, we consider a relatively low mass spectrum of supersymmetric particles. 
Weak scale supersymmetry is one the main classes of new physics signals to be searched for at 
the LHC We adopt several simple SUSY case studies and show that these may cause distortions 
in the Z+jets, TU+jets and tt signal channels. 

In our simulations we adopt the event generator Isajet 7.76 [T] to model Z+jets, PU+jets, 
tt and weak boson pair production at the LHC, along with superparticle production. By 
prescribing a non-zero distribution for W and Z production, Isajet calculates W or Z 
plus one parton emission using exact QCD matrix elements. Further jets are developed via 
the parton shower algorithm. This procedure generates the expected high qriW) and qr{Z) 
distributions. The multi-jet production via this method agrees rather well with overall rates 
for weak boson plus multi-jet production at the Tevatron and is in qualitative accord with the 
merged NLO weak boson production plus parton shower approach advocated in Ref. [8]. 

We use Isajet also for the simulation of signal and background events at the LHC. A toy 
detector simulation is employed with calorimeter cell size Arjx A<p = 0.05x0.05 and —5 < rj < 5. 
The hadronic calorimeter energy resolution is taken to be 80%/ \fE + 3% for \rj\ < 2.6 and 
100%/-\/^ + 5% for the forward region with |?7| > 2.6. The electromagnetic calorimeter energy 
resolution is assumed to be 3%/\^E + 0.5%. We use a UAl-like jet finding algorithm with 
jet cone size R = 0.4 and require that E T (jet) > 50 GeV and \rj(jet)\ < 3.0. Leptons are 
considered isolated if they have Pt{c or /z) > 20 GeV and |?7| < 2.5 with visible activity within 
a cone of AR < 0.2 of TjE^ Us < 5 GeV. The strict lepton isolation criterion helps reduce 
multi-lepton backgrounds from cc and bb production. 

We identify a hadronic cluster with > 50 GeV and \r)(j)\ < 1.5 as a 6-jet if it contains a 
B hadron with Pt(B) > 15 GeV and |?7(-B)| < 3 within a cone of AR < 0.5 about the jet axis. 
We adopt a 6-jet tagging efficiency of 60% and assume that light quark and gluon jets may be 
mis-tagged as 6-jets with a probability 1/150 for E T < 100 GeV, 1/50 for E T > 250 GeV, with 
a linear interpolation for 100 GeV< E T < 250 GeV [9]. 

We will compare expectations for SM calibration reactions against a case study from the 
well-known minimal supergravity (mSUGRA or CMSSM) model, plus two other case studies 
with non-universal soft SUSY breaking terms. In mSUGRA, given a set of parameters 

m , mj/2, A), tan/3, sign(fi), (1) 

one may calculate the entire superparticle mass spectrum, along with sparticle branching frac- 
tions and production cross sections. The Isajet program is used to generate the associated LHC 
collider events. In our work we adopt three case studies listed in Table [TJ The mSUGRA point 
PI is chosen to give a sparticle mass spectrum which yields a modest rate for production of 
real Z bosons from sparticle cascade decays. Point P2 is similar to point PI except that it 
is within the context of non-universal higgs models (NUHM) [10] , wherein the superpotential 
/j, parameter is fixed to a higher value (400 GeV) so that gluino and squark decays to higher 
chargino and neutralino states (which in turn may decay to real Z bosons) are kinematically 
suppressed. Thus, point PI is expected to yield a significant rate for real Z bosons in SUSY 
events, while point P2 is constructed so that Z production in cascade decays is suppressed. 
Finally, point P3 is engineered so as to provide a rather high rate of Z boson production in 
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parameter 


PI 


P2 


P3 


m 


150 


150 


200 




170 


170 


200 


Mi 


170 


170 


120 


^0 


-200 


-200 





tan /? 


20 


20 


20 




262.4 


400 


261.5 


m 5 


432.0 


430.3 


503.8 




416.7 


417.1 


493.4 


HI 7 


256 Q 


302 


352 9 




O04.0 


Q7Q Q 

O ( o.y 


AA1. A 






1 Q9 


944 1 




990 8 


415 2 


294 


Wi 


117 4 


1 26 3 


1 37 6 


m z 4 


9Q0 Q 


414 7 


9Q3 3 


m z 3 


271.1 


404.2 


270.8 




117.7 


126.2 


137.3 




64.4 


66.3 


44.8 


m A 


293.4 


293.4 


327.1 


m h 


108.2 


108.3 


107.7 


a [fb] 


2.2 x 10 5 


2.2 x 10 5 


1.0 x 10 5 



Table 1: Masses and parameters in GeV units for three cases studies points PI, P2 and P3 
using Isajet 7.76 with m t = 171.0 GeV and fi > 0. We also list the total tree level sparticle 
production cross section in fb at the LHC. 



cascade decays. P3 has GUT scale gaugino masses M 2 = M 3 = 200 GeV, while the U(l)y 
gaugino mass Mi is just 120 GeV. The non-universal gaugino masses allow Z 2 — > ZiZ decay, 
and since Z 2 is produced at a high rate in gluino and squark cascade decays, real Z bosons are 
expected to be very common in LHC events for point P3. 

2 Z+jets production at LHC 

To select Z-boson candidate events, we require the following Z Cuts: 

• 2 isolated opposite-sign, same-flavor (OSSF) leptons, 

• 80 GeV < m{U) < 100 GeV, 

• p T (ji) > 100 GeV, 

• Pr{j2, - , in) > 50 GeV, 

where the n jets in each event, ji,j 2 , ■■-,j n , are pt ordered. 
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The staggered jet cuts require the first hard parton emission to be at the highest pt, which 
is modeled by the exact QCD Z + parton matrix element. Further jets at lower p? > 50 GeV 
are generated by the parton shower (PS) algorithm. 

Once the electromagnetic calorimeter and muon chambers are calibrated so that the Z — > 
e + e~ and Z — > peaks are in accord with previous LEP/Tevatron measurements, then 

it will be important to map out the jet activity associated with these events. This will help 
to calibrate the hadronic calorimeter and also to tune the weak boson plus jets Monte Carlo 
programs to the data. Further, the rate for Z — > z/jZ/j+jets production should occur at six times 
the Z —> e + e~+jets rate, so this process gives a good estimate of one of the most important 
backgrounds to the SUSY jets+i?™ lss channel. 

In Fig. [T]we plot the expected jet multiplicity coming from Z — > £ + £~+jets events after Z 
cuts. The red-dashed histogram gives the Isajet prediction for the jet multiplicity from Z+jets 
production alone, while the black-dashed histogram gives the Z plus ti and VV contributions 
(here, V = W or Z). The blue-solid histogram shows the SM contribution plus SUSY from point 
PI. In point PI, the branching fraction of Z 3 -> Z X Z is 11%, Z 3 -> Z 2 Z is 17% and W 2 -> W\Z 
is at 21% , while Z 4 — > ZjZ (i = 1,2) is at 5.7%. However, the squark branching fractions to 
Z 3 , Z 4 and W 2 occur only at the 5-10% level, while gluino branchings into these heavier -ino 
states is less than 1%. Thus, we expect a modest rate for Z+jets production in squark cascade 
decays to Z 3 , Z 4 and H^pJ] for point PI. Point P2 has non-universal Higgs masses, and here 
we have raised the superpotential Higgs mass /x to 400 GeV. This effectively cuts off any squark 
or gluino decays to heavier ino states, so we expect a low rate of Z production in SUSY cascade 
decays. Meanwhile, point P3 has a low GUT scale value of the U(l)y gaugino mass M 1; which 
allows BF(Z 2 — > Z\Z) at 98%. Since Z 2 is produced copiously in gluino and squark cascade 
decays, P3 is expected to yield LHC collider events that are enriched in Z bosons. 

We see from Fig. [T]that the low jet multiplicity events are dominated by Z+jets production 
as expected. As n(jets) increases, ti, WZ and ZZ contributions give elevated predictions of jet 
activity in the "Z+jets" events. In addition, as we step out to higher and higher jet multiplicity, 
the SUSY contribution from point PI becomes ever more important: by n(jets) = 4, the SUSY 
rate is comparable to direct Z+jets production, and by n(jets) = 5, SUSY of PI dominates the 
event rate. Point P3 begins to dominate the Z+jets rate already at n(jets) = 3. Thus, in these 
new physics examples, a direct tuning of SM Monte Carlo programs for Z+jets production to 
data would result in an aberrant tuning, as the data would also include significant portions of 
both SM and new physics events. 

Naively one would expect that the SUSY events could easily be separated from SM back- 
ground by requiring large E™ 1SS - However, in the early stages of LHC running when detectors 
are first being calibrated, the E™ 1SS measurement may not be sufficiently reliable due to a vari- 
ety of detector effects such as mis-calibration, no calibration, under instrumented regions, event 
overlap, beam-gas events, cosmic rays etc.[T2]. In addition, OSSF dilepton events from ti and 
ZZ production will also have E™ 1SS present from neutrinos in the events. 

In Fig. [2J we plot the expected E™ ms distribution arising from Z+jets production (red- 
dashed histogram), Z+jets along with ti and VV production (black-dotted histogram) and all 
SM sources plus SUSY points PI, P2 and P3. The SUSY events typically have much more than 
50 GeV of i?™ lss present. It would then be hoped that, even with initially poor E™ lss resolution, 
the presence of very large E™ 1SS in the Z+jets event sample would alert one to the presence of 
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Figure 1: Jet multiplicity in Z — > ££ events from Z production, tt and VV production and 
sparticle production from SUSY points PI, P2 and P3. 



new physics in this calibration process. 



3 VK+jets production at LHC 

To select W^-boson candidate events, we require the following W Cuts: 

• a single isolated lepton, 

• transverse mass 40 GeV < m T (£, E^ iss ) < 120 GeV, 

• prih) > 100 GeV, 

• p T (j2,--,Jn) > 50 GeV. 

Initially, it will be difficult to establish E™ 1SS as a reliable cut variable, so we omit it from our 
cuts list. However, it is used implicitly in constructing the lepton- missing energy transverse 
mass variable M T (£,Ef ss ). 

One of the first things to check in W — > iug events will be that the transverse mass distribu- 
tion has the requisite shape, including a Jacobian peak at the IV-boson mass. The uncertainty 
associated with the E™ 1SS measurement will lead to broadening of the distribution. As the E™ ss 
resolution is sharpened, the H^-boson Jacobian peak should assume its characteristic form. 
Once confidence is gained that one is really seeing W — > lvt_ events, then the next step will be 
to examine the associated jet activity. The W+jets events will be an important background 
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Figure 2: Missing distribution in Z — > 11 events from Z production, tt and VV (WW, WZ 
and ZZ) production and sparticle production from SUSY points PI, P2 and P3. 



to ti searches, as well as other new physics searches such as SUSY where one looks for single 
lepton plus jets +E^ 1SS events. 

We plot in Fig. [3] the jet multiplicity associated with W —>■ tv^ events after the above cuts. 
The red-dashed histogram shows the contribution from SM IV+jets (without tt contributions) 
events generated from Isajet. The black-dotted histogram includes also the contribution from 
ti and VV events. We also show histograms including PU+jets, ti, VV along with SUSY from 
points PI, P2 and P3. We see that at low n(jets) = 0, 1, indeed most events come from W^+jets 
production. As we move to n(jets) > 2, the events are increasingly dominated by tt production. 
In fact, this also occurs at the Tevatron collider and was proposed as a search strategy for SM 
top events[T3]. The tt contribution to PU+jets events is accentuated at the LHC since the ti 
production rate increases by about a factor of ~ 100 from the Tevatron to the LHC, while W 
and Z production increase by just a factor of 5. We also see that as n(jets) moves out to very 
high multiplicity, n(jets) > 5, it can be the case that SUSY makes a significant contribution 
to the n(jets) distribution. Thus, again, calibrating Monte Carlo W and ti generators to the 
W + n-jets data sample may again lead one askew in that new physics may already be present 
in the actual data event set. 

One way to check early on whether the IV+jets sample at the LHC is in accord with SM 
predictions is to examine the a(W — > eu e )/a(Z — > ee) ratio versus n(jets). In the W/Z ratio, 
the various QCD and PDF uncertainties cancel. In Fig. HI we show the W/Z distribution versus 
n(jets) for the summed SM processes along with SUSY points PI, P2 and P3 at ^/s = 14 TeV. 
From only SM events without ti, the ratio should be at cr(pp —> W) x BE(W — > eu e )/a{pp — > 
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Figure 3: Jet multiplicity in W — > tv^ events from W production, tt and VV production and 
sparticle production from SUSY points PI, P2 and P3. 



Z) x BFiZ — > ee) is around 8.3, as evident in the first bin. This is verified in the first bin: 
n(jets) = 0. In the second bin, n(jets) = 1, the ratio is somewhat lower due to the influence 
of the different M^-cuts versus Z-cuts. At n(jets) > 2, contributions from tt production cause 
an increase in W production relative to Z production, and the ratio steadily climbs with 
jet multiplicity. In the case where Z production is substantial in SUSY events, the SUSY 
contribution eventually will enter, and at high jet multiplicity, the W/Z ratio again starts to 
fall off, indicating the presence of non-SM processes in the calibration sample. This is especially 
so for the Z-rich events from SUSY point P3. 

4 tt production at LHC 

To select tt candidate events, we examine the case where tt — > (Mug) + (bqq'), i.e. one leptonic 
and one hadronic decay. This channel should allow for a simple top mass reconstruction via 
the invariant mass m(jjb) and the cluster transverse mass Mt{M, E^ 1SS ) [T3] , 
We require the following ti 1£ cuts: 

• a single isolated lepton, 

• n(jets) > 4, 

• n(b — jets) = 1, 
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Figure 4: Ratio of W — > e^ e to Z — > ee events at the LHC versus jet multiplicity from SM 
processes (black-solid histogram) and from SM plus SUSY points PI, P2 and P3. 



• transverse mass 40 GeV < m T (£, E™ iss ) < 90 GeV, 

• invariant mass 60 GeV < m(jj) < 90 GeV, and 

• PrCfc) > 40 GeV. 

The jj invariant mass cut selects events where two jets reconstruct a hadronic W. To 
reconstruct the top mass, there is then an ambiguity on which 6-jet is associated with the 
W — > qq' system. Usually we will get this right by choosing the minimum invariant mass 
gained by combining the two jets which make up the W mass with any of the remaining 
jets [15]. 

In Fig. El we plot the invariant mass of the two jets which combine to make up Mw, 
with the remaining jet which yields the minimum invariant mass. The red-dashed histogram 
shows the resulting distribution from ti production, and a clear peak is seen just below the 
value of m t . The high invariant mass tail comes from a wrong assignment of the additional 
jet to the m(jj) = My/ cluster (here the correct parton is usually too soft to satisfy the jet 
finding criteria, or ends up overlapping with another jet). We also show the ti plus the l^+jets 
contribution as the black-dotted histogram, which essentially overlaps with the ti plot since 
the H^+jets background is very low for this set of cuts. We also show in Fig. El the invariant 
mass construction when the SUSY contributions from point P2 are added in. Here, we again 
see the distinct top peak, but now the high m(jjb) shoulder is much higher. In this case, the 
reconstructed top mass distribution would suffer significant shape distortion on the high end. 
A tuning of the ti Monte Carlo to LHC top quark data could again result in a mis-calibration. 
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The shape distortion can be simply characterized by the peak height to the shoulder height. In 
the SM, the ratio of heights is ~ 2, while adding in a SUSY contribution gives a ratio of ~ 1.5. 

In Fig. [6] we show the cluster transverse mass Mt(M, EJp lss ). This distribution is formed 
from the isolated lepton combined with the jet which yields the smallest invariant mass and then 
constructing the cluster transverse mass with the E™ 1SS variable. The distribution in the ideal 
case would be a continuum bounded from above by m t . Here, we see that the ti contribution has 
a tail extending to high cluster transverse masses owing to mis-combination of jets with leptons, 
and also from additional sources of E™ 1SS , such as neutrinos and jet energy mis-measurement. 
The distribution from ti plus SUSY with point P2 has some widening associated with it, and 
so could again lead to Monte Carlo mis-tuning. 



tt, single lepton cuts 




M T (bl,E T )(GeV) 

Figure 6: Distribution in cluster transverse mass m<r(W, E^ 1SS ) from ti and IU+jets production, 
and from SUSY production with m = 150 GeV, m 1/2 = 170 GeV, A = -200, tan/3 = 20, 
/i > and m t = 171 GeV. 

One of the main ways to distinguish if the ti signal is really from ti is to compare results from 
the \i channel to results from the dilepton channel, where both t and i decay semi-leptonically. 
We will invoke the following ti It cuts, taken from Ref. [T6] : 

• two isolated leptons, 

• n(jets) > 2, 

• n{b — jets) = 2, 
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Table 2: Cross sections in fb from ti cuts for 1£ events and 2£ events from SM sources and 
for SM plus contributions from mSUGRA points PI, P2 or P3, along with the single lepton to 
dilepton ratio. 
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0.014 
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Table 3: Cross sections in fb for OS and SS dilepton events from SM sources and for SM plus 
points PI, P2 or P3, along with the SS to OS dilepton ratio. 



• p T (£ 2 ) > 20 GeV, and 

• Ef ss > 40 GeV. 

For now, we will not distinguish whether or not the dileptons are opposite-sign (OS) or same- 
sign (SS). The ratio of single-to-dilepton events from ti and other SM sources should be a very 
predictable quantity. However, in SUSY events, there is a high propensity to generate states 
with large isolated lepton multiplicity due to the gluino and squark cascade decays. Thus, if 
there is SUSY contamination of the ti signal, we would expect enhanced dilepton signal rates 
compared to single lepton rates. We illustrate this in Table [2, where the single lepton and 
dilepton rates are listed in fb after tt-cuts for single (C3 cut) and dilepton (C4 cut) events. We 
see a decrease in the l£/2£ ratios due to SUSY contamination. 

Another possibility is to separately consider OS and SS dileptons. SS dileptons are expected 
to be rare at LHC, although they can arise from ti production where t — > b£ + vi and i — > bqq', 
followed by b — > c£ + vi decay, where the lepton from b decay is isolated, and from B° — B° 
mixing. Further sources come from secondary processes such as W ± W ± production, which we 
do not simulate here. However, in SUSY events, SS dileptons can arise from gg production, 
where each gluino decays to W~i + qq' for instance [Tf] . The SS dileptons can also easily arise 
from gq and qq production. In Tabled we list OS and SS dilepton rates in fb from SM sources 
along with SUSY points PI and P2. Here, we see that the SS/OS ratio is greatly enhanced 
when SUSY events are present. Further, since LHC is a pp collider, we expect an asymmetry 
in SS dilepton production when squarks are produced at large rates. This is also seen in Table 
[31 where a measurable ++ to asymmetry is observed due to SUSY contamination. 
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5 Conclusions 



One of the first orders of business for the LHC experiments will be to "re-discover" the Standard 
Model, i.e. to make sure that familiar processes like QCD dijet production, lU+jets production, 
Z+jets production and ti production all occur as expected. In fact, these processes will likely 
be used as calibration processes, to calibrate the detectors and to tune Standard Model Monte 
Carlo programs. In this paper, we have emphasized that low mass SUSY may already make 
significant contributions to these familiar SM processes, these SUSY contributions, if not taken 
into account, could cause detector mis-calibration, or mis-tuning of SM Monte Carlo codes. 

It will be imperative to check that the SM calibration reactions have all of the properties 
expected from these reactions. We point out several ways in which Z+jets, VU+jets and ti 
production can be checked for contamination from non-SM processes like low mass SUSY. 

• In Z+jets events, even though initially E™ ss will not be well-measured, there should still 
be an obvious excess of E™ ss from SUSY contamination. 

• In IU+jets events, we advocate plotting the W/Z ratio versus jet multiplicity. The SM 
prediction should grow with jet multiplicity due to the presence of ti contributions. How- 
ever, SUSY contamination is likely to upset the expected ratios, especially at high jet 
multiplicity. 

• In ti events, we advocate comparing the shape of the trijet distribution which reconstructs 
m t against SM expectations. In addition, the dilepton-to-single-lepton event ratio can be 
sensitive to SUSY contamination, as can the presence of a high rate of same-sign dileptons. 

Only if all the properties from the SM processes accurately match the data can they be 
used with a high degree of confidence for calibrations. 
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